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Abstract. We have investigated the turbulent mean-field dynamo action in protoneutron stars that are subject to convective and 
neutron finger instabilities. While the first one develops mostly in the inner regions of the star, the second one is favoured in 
the outer regions, where the Rossby number is much smaller and a mean-field dynamo action is more efficient. By solving the 
mean-field induction equation we have computed the critical spin period below which no dynamo action is possible and found 
it to be ~ Is for a wide range of stellar models. Because this critical period is substantially longer than the characteristic spin 
period of very young pulsars, we expect that a mean-field dynamo will be effective for most protoneutron stars. 
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1. Introduction 

The origin of the strong magnetic fields in neutron stars is still 
a matter of controversy. The magnetic fields inferred from the 
pulsar spin-down data ranges from ~ 5 x lO^'^ to lO^G 
but these values are representative of the global magnetic con- 
figuration rather than of the fine magnetic structure near the 
stellar surface. X-ray spectra of some pulsars have recently 
started to provide a closer look at the magnetic field strength 
near the neutron star surface and the absorption features in the 
spectrum of IE 1207.4-5209, for instance, have been used to 
estimate a strong surface magnetic field, Bs ^ 1.5 x 10^*G 
(Sanwal et al. 2002). This is to be contrasted with the dipo- 
lar magnetic field estimated from the spin-down rate of this 
pulsar: B<i - (2 - 4) x IQi^G (Pavlov et al. 2002), a value 
rather typical for a radio pulsar of ^ 0.2 — 1.6 Myr. As another 
example, Becker et al. (2002) have reported the presence of 
an emission line in the X-ray spectrum of PSR B 182 1-24 that 
could be interpreted as cyclotron emission from a corona above 
the pulsar's polar cap. The line would be formed in a magnetic 
field Bs ^ 3 X 10^^ G, approximately two orders of magnitude 
stronger than the dipolar magnetic field inferred from the spin 
evolution. Both measurements provide evidence that the local 
magnetic fields at the neutron star surface can be well above the 
dipolar one responsible for the secular spin-down of pulsars. 

Observations of radio emitting pulsars also exhibit a 
distinction between the dipole and surface magnetic fields. 
Recently, Gil & Mitra (2001) and Gil & Melikidze (2002) have 
argued that the formation of a vacuum gap in radio pulsars is 



possible if the actual surface magnetic field near the polar cap 
is very strong, B^ ^ lO^^G, irrespective of the magnetic field 
measured from the spin evolution. Furthermore, the presence 
of a strong magnetic field with a small curvature (< 10*' cm) 
can account for the radio emission of many radiopulsars that 
lie in the pulsar graveyard and should be radio silent (Gil & 
Mitra 2001). This growing number of evidences for a complex 
structure of the magnetic field at the surface of neutron stars 
suggests that this may represent a general property of pulsars. 

Magnetic fields with different strengths on different length- 
scales can be explained naturally if they are generated through 
a dynamo mechanism driven by turbulent motions. Indeed, it 
is generally accepted that protoneutron stars (PNSs) are sub- 
ject, shortly after their birth, to hydrodynamic instabilities in- 
volving convective motions (Epstein 1979, Livio et al. 1980, 
Burrows & Lattimer 1986) and that these can last ^30 — 40 
s (Miralles et al. 2000, 2002). We here show that, under suit- 
able conditions, turbulent motions can generate magnetic field 
via dynamo action and that a mean-field dynamo can be opera- 
tive together with small-scale dynamo processes (Thompson & 
Duncan 1993; Xu & Busse 2001). 

2. Convection in PNSs 

Hydrodynamic instabilities in PNSs can be driven by either lep- 
ton gradients (Epstein 1979) leading to the so-called "neutron- 
finger instability" (Bruenn & Dineva 1996), or by negative 
entropy gradients which commonly observed in simulations 
of supernova explosions (Bruenn & Mezzacappa 1994, 1995; 
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Rampp & Janka 2000) and evolutionary models of PNSs (Keil 
& Janka 1995; Keil et al. 1996; Pons et al. 1999). This latter 
instability is usually referred to as the "convective instability", 
although both instabilities involve convective motions. 

It has been calculated that these instabilities will first de- 
velop in the outer layers containing ^ 30% of the stellar 
mass, with the convectively unstable region surrounded by the 
neutron-finger unstable region, the latter involving therefore a 
larger portion of the stellar material. After a few seconds the 
two unstable regions move towards the inner parts of the star 
and after ~ 10 s from the initial development, more than 90% 
in mass of the star is hydrodynamically unstable. At this stage 
the stellar core has become convectively unstable but it is still 
surrounded by an extended neutron-finger unstable region. In 
the ~ 20 s that follow, the temperature and lepton gradients are 
progressively reduced and the two unstable regions begin to 
shrink, leaving the outer regions of the star. After ~ 30 s, most 
of the PNS is stable and the instabilities disappear completely 
after ~ 40 s (Miralles et al. 2000). 

During this period the PNS is opaque to neutrinos and the 
turbulent mean velocity can be estimated within the mixing- 
length approximation. The largest unstable length-scale is then 
of the order of the pressure length-scale, L = p\dp/dr\~^, and 
the effective flow velocity in this scale, , can be estimated as 

K, L/t^, where is the growth-time of instability which is 
of the order of the turnover time in the scale L (Schwarzschild 
1958). In some cases, can be computed by equating p?;^, 
with p the rest-mass density, to the total neutrino energy flux 
(Thompson & Duncan 1993). While the two approaches yield 
similar results at the peak of the convective instability, the first 
one is expected to be more accurate especially when not all of 
the energy is transported by turbulence (e.g. when the insta- 
bility is not fully developed), when the temperature gradient is 
not significantly super-adiabatic, or when the lepton gradient is 
small. 

In the convectively unstable region, the growth-time of in- 
stability is 



t2 3^^ T 



(1) 



where Tc is the growth-time of convection, g is the gravitational 
acceleration, A VT is the difference between the actual and the 
adiabatic temperature gradient, and j3 is the coefficient of ther- 
mal expansion. Except during the last stages of the unstable 
phase, when the entropy gradients have been washed out, the 
convective instability grows on a short dynamical timescale. 
Miralles et al. (2000), have estimated this to be Tc ^ 0.1 — 1 
ms, from which we derive the a mean turbulent velocity in the 
convectively unstable region ti^ ~ 10^ — 10^ cms~^. 

In the neutron-finger unstable region, on the other hand, the 
lepton number gradients dominate over the temperature gradi- 
ents and we can estimate the growth-time as 

4 ~ 4- ~ \9^wn > (2) 

where 6 is the coefficient of chemical expansion, and Y = 
{ue + n„)/n is the lepton fraction with rig, n^, and n be- 
ing the number density of electrons, neutrinos, and baryons. 



respectively. The estimated growth-time in the neutron-finger 
unstable region is a couple of orders of magnitude longer than 
the one for the convective instability, i.e. t„/ 30 — 100 ms 
(Miralles et al. 2000), thus yielding a mean turbulent velocity 
t;^ ~ (1-3) X lO^cms-^ 

The existence of two unstable regions with substantially 
different mean velocities is the most significant difference with 
the model by Thompson & Duncan (1993) in which the whole 
PNS was assumed to be convectively unstable with turbulent 
velocity ~ 10* cm s~^. Furthermore, a longer growth-time 
(and hence turnover time) is what promotes the efficiency of 
the mean-field dynamo in the neutron-finger unstable region. 
We recall, in fact, that PNSs are likely to rotate and although 
the initial spin rates of pulsars are not well constrained by ob- 
servations, they are believed to be around 100 ms (Narayan 
1987). As a result, the Rossby number, Ro = P/t^, with P 
being the PNS spin period, can take substantially different val- 
ues in the two unstable regions. In particular, Ro ^ 100 in 
the convectively unstable region and the influence of rotation 
on the turbulence is therefore weak. As a consequence, and as 
already pointed out by Thompson & Duncan (1993), the mean- 
field dynamo will not operate efficiently here. In the more ex- 
ternal regions unstable to neutron-fingers, on the other hand, 
Ro ^ 1, turbulence can be strongly modified by rotation and 
this favours the efficiency of a mean-field dynamo. Of course, 
in both regions turbulent magnetic fields can also be generated 
by small-scale dynamo driven by turbulent motions. 

To investigate more quantitatively the efficiency of a mean- 
field dynamo action, we have modelled the PNS as a sphere 
of radius i? having two spherical turbulent zones with substan- 
tially different properties and separated at Rc. The inner parts 
(r < Rc) correspond to the convectively unstable region, while 
the outer ones (Rc < r < R)io the neutron-finger unstable re- 
gion. The boundary between the two regions moves inward on 
a timescale comparable to the cooling timescale (i.e. 1 — 10 
s), much longer than the turnover time for both instabilities. 
Hence Rc is a slowly varying variable and it is sufficient to 
consider a sequence of PNS models differing only for Rc. 

3. Dynamo action in PNSs 

Numerical simulations indicate that the turbulence in PNSs 
will be non-stationary, developing rapidly soon after the col- 
lapse, reaching a quasi-stationary regime after a few seconds, 
and then progressively disappearing as the temperature and lep- 
ton gradients are removed. Because the characteristic cooling 
timescale for the PNS exceeds both t„/ and Tc, the turbulence 
can be treated adiabatically (this assumption will cease to be 
accurate as the instabilities are progressively suppressed). In 
this case, the mean-field induction equation for a turbulent, 
magnetised and conducting plasma can be written as 



— = yx{vxB + aB) 



V X (r?V X B) , 



(3) 



where is the turbulent magnetic diffusivity, a is a pseudo- 
scalar measuring the efficiency of the dynamo (the "a- 
parameter"). Here, v is the velocity the ordered fluid motion, 
which we assume to follow a simple law, v = x r, but 
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allow for the differential rotation often observed in numeri- 
cal simulations (Zwerger & Miiller 1997; Rampp et al. 1998; 
Dimmelmeier et al. 2001). As customary in dynamo theory, we 
express this differential rotation in terms of the radial coordi- 
nate r and through a simple quadratic law, i.e. 

n{r) = n„ + r^n, , (4) 

where the coefficients il^ = il{r 0) > and il^ are not nec- 
essarily chosen so as to satisfy the Rayleigh stability criterion. 
Boundary conditions for the magnetic field need to be speci- 
fied at the stellar surface, where we impose vacuum boundary 
conditions, and at the centre of the star, where we impose the 
vanishing of the toroidal magnetic field. 

As discussed in Section 2, all of the PNS undergoes tur- 
bulent motions but with properties that are different in the in- 
ner parts (0 < r < Rc), where fast convection operates, from 
those in the outer parts (Rc ^ r < R), where the neutron- 
finger instability operates. To model this in a simple way, we 
assume the relevant physical properties of the two regions to 
vary in a smooth way mostly across a thin layer of thickness 
AR = 0.025R. More precisely, we express 77 as 

V^Vc + iVnf - Vc) {1 + erfiir - Rc)/AR]}/2 , (5) 

where 77^ and ri„ f are respectively the turbulent magnetic diffu- 
sivities caused by the convective and neutron-fingers instabil- 
ities, and erf is the "error function". Using (|5}, r; w 77^ in the 
convectively unstable zone (Rc — r AR), while 77 w r/„/ in 
the neutron-finger unstable zone (r~Rc 3> AR). Furthermore, 
because in typical PNSs 77„/ <C 77c, we have here chosen 
rjnf/rjc ~ 0.1. Similarly, we have modelled the a-parameter 
as being negligibly small in the convectively unstable region 
and equal to q;„/ in the neutron-finger unstable region, i.e. 

a{r,B) =a^f cose {\ + erf[{r ~ Rc)/ AR]} /2, (6) 

where the angular dependence chosen in (|6|l is the simplest 
guaranteeing antisymmetry across the equator. 

We recall that in a rotating turbulence with lengthscale £ 
and moderate Rossby number, a « — 17^^V ln(pi;^ ) (Riidiger 

6 Kitchatinov 1993). In PNSs, however, the pressure is deter- 
mined by degenerate neutrons and can, in a first approximation, 
be expressed with a simple relation of the type p oc , where 

7 = 5/3 for a non-relativistic neutron gas and 7 = 4/3 for a 
relativistic one. As a result, the density lengthscale is compa- 
rable to the pressure one L and, as mentioned in Section 2, to 
the lengthscale of the instabilities. This introduces a great sim- 
plification since we can express the isotropic turbulence in the 
neutron-finger unstable zone simply as Q!„/ « Q.L. 

4. Numerical results 

The induction equation (|3jl with 77 and a given by (|5} and 
(|6jl has been solved with a numerical code employing finite- 
difference techniques for the radial dependence and a polyno- 
mial expansion for the angular dependence (see Bonanno et al. 
2002 for details). The simulations reported here use 30 spheri- 
cal harmonics and about 40 grid points in the radial direction; 
details on the numerical procedure will appear elsewhere. 
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Fig. 1. Critical period as a function of the differential rotation 
parameter The two curves refer to different values of i?c, with 
the solid parts corresponding to a -dynamo and the dashed 
parts to a af2-dynamo. 



Assuming Vj^Lj?, = 10^^ cm^ s^^ and progressively vary- 
ing the velocity field (0} in terms of the differential rotation 
parameter q = R^il,-^/V,{r = R), we have solved equation (|3} 
to determine the critical value corresponding to the marginal 
stability of the dynamo. Given a certain amount of differential 
rotation, in fact, the seed magnetic field will grow if a„/ > 
and instead decay if «„/ < . The different types of dynamo 
can be distinguished according to whether the differential rota- 
tion is small and the evolution of the magnetic field stationary 
(i.e. a a^-dynamo), or viceversa (i.e. a ail-dynamo). 

Since a„/ ~ ilL, the critical value oq effectively selects 
a critical value for the spin period, = 2ttL /uq, such that 
magnetic field generation via a mean-field dynamo action will 
possible only if the stellar spin period is shorter than the critical 
one. In Fig. 1 , we plot the critical period as a function of the dif- 
ferential rotation parameter q. Note that q < and q > corre- 
spond to situations in which the stellar surface rotates faster and 
slower than the centre, respectively (values g < — 1 correspond 
to a counter-rotation and may be not physically relevant). 

As shown in Fig. 1, a stationary a^-dynamo dominates the 
magnetic field generation process for \q\ < 1, while a afl- 
dynamo is more efficient for \q\ > 1. This latter case will be 
characterised by a magnetic field of oscillating strength but, 
given the large differential rotation required, it may be diffi- 
cult to achieve in practice. Hence, the a^-dynamo appears to 
be the most likely source of magnetic field generation via dy- 
namo processes in PNSs. 

The critical spin found here is in general rather long and 
for a PNS rotating uniformly (i.e. q = 0),a mean-field dynamo 
will develop if P < 1 s when Rc/R = 0.6 and if P < 1.5 
s when Rc/R = 0.3. This difference is due to the fact that a 
PNS with a more extended neutron-finger unstable region can 
rotate proportionally more slowly while maintaining the same 
dynamo action. Furthermore, if the star rotates differentially 
with \q\ w 1, Pg is further reduced, being approximately 20% 
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Fig. 2. Toroidal (B^) and the poloidal (Bp) magnetic field fines. 
Solid and dashed contours correspond to positive and negative 
values, respectively. The dot-dashed line marks the position of 
Rc = 0.6R. 

shorter. As a result, only PNSs with P > 1 — 1.5 s will not de- 
velop a turbulent mean-field dynamo. Such slow rotation rates 
should be rather difficult to achieve if angular momentum is 
conserved during the collapse to a PNS and, indeed, observa- 
tions suggest that the initial periods of pulsars are considerably 
shorter than the critical period obtained here. We expect, 
therefore, that a turbulent mean-field dynamo will be effective 
during the initial stages of the life of most PNSs. 

In Fig. 2 we show the toroidal (-B,/,) and poloidal {Bp) mag- 
netic fields of a typical PNS model. Note that both components 
are generated in the outer neutron-finger unstable region, but 
turbulent diffusion produces a magnetic field also in the in- 
ner convectively unstable region, although this is considerably 
weaker. The toroidal magnetic field tends to concentrate near 
the polar regions whereas the poloidal one is more evenly dis- 
tributed in latitude. The numerical calculations also show that 
if Igl < 1 and the field is generated by the a'^ -dynamo then 
B^/Bp - 10, while B^/Bp ~ 100 - 200 if \q\ > 1 and the 
afi-dynamo generates the magnetic field. Both results suggest 
that the internal magnetic fields in neutron stars could be sub- 
stantially stronger than the observable surface ones. 

5. Conclusions 

We have calculated the turbulent mean-field dynamo during the 
turbulent instabilities that are expected to accompany the early 
stages of the life of a PNS. For - 30 - 40 s, in fact, the PNS is 
subject to two substantially different instabilities, with a con- 
vective instability active in the inner regions of the star and 
a neutron-finger instability being more efficient in the outer re- 
gions. The turbulent motions are more rapid in the convectively 
unstable zone, where the Rossby number is large and the a- 
parameter is likely not to be significant. In the neutron-finger 
unstable region, on the other hand, the turnover time is consid- 
erably longer, the Rossby number small, and the a-parameter 
is sufficiently large that a mean-field dynamo can operate. 

The occurrence of a dynamo depends sensitively on the 
stellar rotation rate since this can influence the turbulent mo- 



tions if sufficiently high. Our simulations show that even rela- 
tively slowly rotating PNSs can be subject to a dynamo action 
in the neutron-finger unstable region, with the -dynamo be- 
ing the most efficient mechanism. The critical value of the spin 
period below which the dynamo is suppressed has been found 
to be ~ Is for a wide range of models and is essentially 
larger than the characteristic spin period of young pulsars as 
inferred from observations. As a result, a turbulent mean-field 
dynamo can be effective in the early stages of the life of most 
PNSs. The magnetic field produced in this way is concentrated 
mostly where the generation occurs, but turbulent diffusivity 
transports part of it also in the inner regions of the star 

We have here considered the generation of axisymmetric 
magnetic fields, but the -dynamo could also lead to the gen- 
eration of non-axisymmetric magnetic fields, with the critical 
spin period being not very different from the one discussed in 
the axisymmetric case (Riidiger et al. 2003). The generation of 
non-axisymmetric magnetic fields and the effect of saturation 
will be considered in a forthcoming paper 

We note that the production of the mean magnetic field by 
a turbulent dynamo is accompanied by the generation of small- 
scale magnetic fields which are likely to be stronger than the 
mean one and that have not been considered here. However, 
nonlinear effects associated to these small-scale fields can in- 
fluence the dynamo and need to be taken properly into account. 
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